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COMPUTER PROGRAM FOR SOLVING NINE -GROUP DIFFUSION 
EQUATIONS FOR CYLINDRICAL REACTORS 
By James W. Miser , Robert E. Hyland, and Daniel Fieno 


SUMMARY 

This report presents a method for determining the critical size of 
a cylindrical reactor by a one-dimensional group-diffusion solution in 
h the radial plane. Hie effect of leakage in the axial direction is taken 

into account by prescribing values of axial leakage based on assumed 
flux levels. The method is based on dividing the neutron energy spectrum 
into nine groups and the reactor into four concentric cylinders (or re- 
gions) . The nuclear cross sections used for each group are the average 
values in each energy group. 

A computing machine program for an IBM 650 computer (plus attach- 
ments specified herein), which performs the solution of the diffusion 
equations given, and a method for using the program for less than nine 
groups and four, or less, regions are included. 

INTRODUCTION 

The subject program for the IBM type 650 magnetic drum data-processing 
machine provides a method for determining the critical size of a cylindri- 
cal reactor on the basis of age-diffusion theory. The age-diffusion equa- 
tion is approximated by group-diffusion equations wherein the entire 
neutron energy spectrum is divided into nine energy groups selected on 
the basis of neutron cross-section characteristics* The neutron cross 
sections and diffusion coefficients for each energy group are obtained 
by a suitable averaging procedure performed prior to the subject program. 
The program is written in such a way that fissions can occur in any, or 
all, of the energy groups and the neutrons may be born into any, or all, 
of the energy groups . 

The reactor consists of four concentric cylinders, herein called re- 
o i on . . Any, or all, of the inner three regions may contain fissionable 
material, but the outer region may not. In the axial direction, each 
region is assumed, to be uniform in composition. 
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By proper adjustments of the input information, the number of energy 
groups and/or regions can be reduced. 

The subject program was written in the language of IBM SOAP II 
(Symbolic Optimal Assembly Program) and then converted by the IBM 650 
computer into basic machine language. The machine time required for 
the solution of one case is approximately 6 minutes. Thus, in approxi- 
mately 18 minutes three values of the effective multiplication factor 
can be plotted as a function of one geometric variable to determine the ^ 

approximate critical size of the reactor in question. 

This report presents the method of analysis and a discussion of the 
calculating procedure used to determine the critical size of a reactor. 
Information is included on the use of the machine program given and the 
use of the program for less than nine energy groups for four, or less, 
regions, as well as an outline and tabulation of the machine program. 

The method of obtaining the necessary input data is similar to that of 
reference 1, but substitutions of other methods for obtaining the re- 
quired input can be made at the discretion of the reader. References 2 
to 5 may be of use in understanding the programming and operating pro- 
cedures used with this program. 


METHOD OF ANALYSIS 

The program calculates the effective multiplication factor k e ff 
of a reactor by means of the equations and procedure given herein. On 
the basis of three or four such solutions for various values of one of 
the geometric or physical variables, the critical size of the reactor 
can be obtained. 

The solution for kg^-p is begun by first assuming a radial distri- 
bution of neutrons born and then dividing the neutrons into energy 
groups on the basis of the division of the fission spectrum. From this 
assumed radial and energywise flux distribution of neutrons born, a new 
radial and energywise flux distribution can be calculated. On the basis 
of the new flux distribution, a calculation of the neutrons born as a 
result of the fissions in each energy group is made. The reactor power 
is proportional to the number of neutrons born. The new value of reactor 
power is divided by the previous one to determine a corresponding value 
of the effective multiplication factor k eff . The new radial distribu- 
tion of neutrons born then replaces the originally assumed distribution, 
and this cycle of calculations is continued until the value of kg.^ is 

stabilized; or, as considered in the subject program, Ak e ff approaches 
zero. 
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The equations used to make these calculations are developed in the 
following sequence: 

(1) The steady- state group- diffusion equation is written for each 
energy group employing group-averaged cross sections, 

(2) The group- diffusion equation is then written in finite differ- 
ence form for selected radial points in the reactor, 

(3) A procedure for solving equations at each radial point is de- 
scribed for determining the neutron flux at each point for the energy 
group in question. 

(4) For the radial flux distribution just calculated, the equation 
is given for computing the radial distribution of neutrons born out of 
a particular energy group, 

(5) Using the summation of neutrons born at each radial point, 
equations are given for calculating the reactor power and k e ff° 

The value of the results obtained from the group-diffusion calcula- 
tions may in some cases be improved with an increase in the number of 
energy groups. For this reason, the number of energy groups selected is 
nine, which is the largest number of groups that could be accommodated 
by the computing machine for which the program was designed. These nine 
groups are numbered a to i, and any one group generally will be re- 
ferred to as the n th group. Group a designates the group of highest 
neutron energy, and group i designates the group of neutrons at thermal 
energy. 

At any point within the reactor, the steady- state equation can be 
written in the form: 


Neutron loss = Neutron production (l) 

When the case of a unit volume within the reactor and neutrons of a 
particular energy group n is considered, the terms of the equation 
can be subdivided as follows (all symbols are defined in appendix A) : 

Neutron loss per unit time : 

2 

(1) Diffusion of neutrons from volume element = -DqV <p n 

(2) Neutron absorption = Za,n<Pn 

(3) Slowing out of group n = 2q, n <Pn 

(4) Neutron leakage from reactor in axial direction = DnB^ri 
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Neutron production per unit time : 

(l) Fission neutrons born into group n = a n pt where 
i 

b, = v y* <p 

M t / j f^irn 


(2) Neutrons slowing down into group n = 2^ 

If the previous terms are substituted in equation (l), the general 
diffusion equation becomes 


-D n V^ n + ( Z a ,n + Z q,n + V^n = + 2, 


'q,n-l < ’ > n-l 


( 2 ) 


Grouping the terms within the parentheses into one term, 2 n <p n reduces 
equation (2) to the form 


-V % + Vn = °n*t + Z ~ 


q,n-l n-1 


( 3 ) 


For a given unit volume in the reactor, the diffusion equations for the 
nine energy groups can be written 


(-D A + 2cp) a = a & p t + 0 
(-DV 2 <p + Z<p) b = +2 q ^ a (p a 

(-DV2(p + Zcp) c = a c p t + 2^ b c^ L 


(-DV 2 cp + 2<p) i = 0 + 2^ h <p h 


( 4 ) 




For the cylindrical geometry selected, the flux cp is assumed to 
be a function of the radius r only; therefore, 


o d^tp 1 d<P 

V^cp = — - + — 


dr 


r dr 


(5) 


The terms on the right side of equation (5) may be approximated by ex- 
pressing them in terms of finite differences based on values at adjacent 
points within a region. For this reason, each region is subdivided into 
cylinders of equal thickness h except for the cylinder at the center 
of the reactor, the two cylinders on either side of an interface, and 
the outermost cylinder (see fig. l) . Figure 1 presents the assumed 
configuration of the reactor and indicates the location of the numbered 
points. The radial dimensions of the four regions (i, II, III, and IV) 
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may vary independently; in other words , the distance h between two ad- 
jacent radial points in a region may vary from region to region. Thus, 
the value of h is specified by means of the general subscript m e 

At a point k the derivatives of equation (5) can be written in 
terms of finite differences as 


/cW>\ k \+l ~ \-l 

V4rj k 2h m 


( 6 ) 


/dM = ^ k+1 ~ 2 ^k + ^k-l 

W 2 /k h m 

Substituting equations (6) and (7) into equation (5) and simplifying 
result in the following equation: 


V 2 cp 


E <P 
k-lk-1 


2r k\ 


- Vk+1 


r, h 2 
k m 


( 8 ) 


where 


E. 


'k-l 


= r k “ 


lx m 

2 


E k = 


+ 


2 


Substituting equation (8) into equation (3) and grouping the coefficients 
of <P k give 


D 

n,m 

r k^ 

= “n^tjk + 2 q J( n-l,m <,> k,n-l 
Further simplification gives 

Ek-A-l^n + r k^n ; m < * ) k,n + ^k+l^n = W k,n 


^k-l^k-l.n + 


•2r> - 


r,h 2 S \ 

k m n,m \ 


D 


<P 


k,n 


E k^k+l,n 


n ,m 


( 9 ) 
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where 


h 2 Z 

jnn^ 

n D n,m 

r k h 2 
k m 


W k,n " D^ m (“n^k + 2 q,n-l,m<Pk,n-l) 


Equation (9) for the n th group is then of the form 


w 


C l,k*k-1 + C 2 ^ k <p k + C 3A cp k+1 = w k (10) 

Equation (9) is applicable at all points in the reactor except where 
special conditions are imposed, such as points 1, 11, 21, 31, and 38 in 
figure 1. Points 1 and 38 can he treated as special cases of equation 
(9)$ however, points 11, 21, and 31 at the interfaces of adjacent regions 
cannot be solved by equation (9) because the value of G n is dependent 

on the cross sections of a particular region. 

In the case of point 1, which is a distance of hj/2 from the 
center of the reactor, it is assumed that <p = <p Q ; therefore, by sub- 
stituting hj/2 for h m and for <Pq, equation (9) becomes 



" 8 D^ t ^^,1 + ^,11-1,^1,11-1) 


(The subscript I refers to the central region (fig. 1(a)).) Dividing 
this equation by hj/4 and simplifying gives 


-3 - 


h T 2 T 

I n.I 


2 + 5<p 2,n “ " 2D n ^j (“fc^t,! + Z q,n-l,fPl,n-l) 


( 11 ) 


which 


, in the form of equation (10) for the n th group, is 


C 2,l (p l + C 3,l <(> 2 = W 1 


£TT“ 
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For the solution of equation (9) at point 38, it is assumed that 
T 39 = 0 and that 0t,38 = 0; therefore, equation (9) can be written as 
follows : 


^37^37, n + 


r 38 u 38,n v 38,n 


,2 

' 38 IV 


D 


n,IV 


^q,n-l,IV < P38,n-l^ ( 12 ) 


or , after the general form of equation ( 10 ) for the group, 

G l, 38^37 + C 2, 38^38 = w 38 

In order to determine an equation that is suitable for computing 
the flux at the interface of two regions (i.e., points 11, 21, and 3l), 
the requirement that the net neutron current (-D grad <p) be the same 
for each region at the interface is considered. At point 11, for 
example, 


-D t (grad <p). 
I,n^° 


, n , 11 


D II,n (grad <P) 


II,n,ll 


This can be expressed as 


-D (~) 

I,n\dr/ T 
' I,n,ll 


-D. 


II 


,n\dr/ 


II,n,ll 


(13) 


The value of (d<p/dr)j n ^ can be determined approximately by the fol- 
lowing equation, representing the first two terms in a Taylor's series 
expansion of d<p/dr at point 10 : 



I,n,ll 


* 



+ h. 


kdr 2 /n ,10 


Substituting for the derivatives their equivalents as written in equa- 
tions ( 6 ) and (7) gives 


q<p\ 


< P 


11 , n 


9 ,n 


2hi 


+ h T 


^11^ 


2<P 


10^11 


hi 



which, simplified, is 


/d«P\ 

W) 


I,n,ll 


5 b.l,n - 4 (p 10 ,n + »9,n 
2 h T 


( 14 ) 
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By a similar procedure, 




, dr 7ll,n,ll 


3® - 4<P + <P 

11, n 12, n 13, n 


2h 


II 


( 15 ) 


Rearranging equation (13) in the form 

(M) _ = o 

^' 1 , 11,11 D I,n ^11^,11 
and substituting equations (14) and (15) give 


3T _ 44> + <p 

11, n " ' 


2hj 


l>I,n 


3<l> 


y xx 10, n 9,n II, n I 11, n 12, n 15, n 


-4 <P + <P 

12,r - 1 

2hn 


Then, simplifying this equation gives 


<p Q - 4<p.._ + 3(1 + a T )<P,, - 4a <P + a T <P, „ = 0 

9,n 10, n I,n 11, n I,n 12, n I,n 13, n 


(16) 


where 


l I,n 


h T D TT 

I II, n 
h II D I,n 


At each interface (or k = 11, 21, 3l), equation (16) can be written in 
general form for the group as follows : 


•k-2 + d i,kVi + 


d 2,k\ 


d 3,k\+l + 


d. . = 0 

4,k k+2 


(17) 


As a result of this development, equations (9), (ll), (12), and 
(16) make it possible to write one equation at each of the 38 radial 
points in the reactor in terms of the flux with coefficients based on 
the nuclear cross-section data for each energy group. The flux coeffi- 
cients on the left side of each equation for the n^h group form a 
matrix (designated in the machine program as the left matrix) in which 
each row represents an equation at a particular radial point and each 
column represents the flux at a particular point as indicated: 
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9 


3 

H 

-I 

1 


VJ 

I 


C 2,l G 3,l 
C l ,2 C 2,2 C 3,2 


C l,3 C 2,3 C 3,3 


C l,10 

C 2,10 

C 3,10 



1 

d l,ll 

d 2jll 

d 3,ll 

d 4,ll 



C l,12 

C 2,12 

C 3,12 


C l,36 

C 2,36 

C 3,36 



C l,37 

C 2,37 

C 3,37 



G l,38 

C 2,38 


The two simultaneous equations represented by the first two rows 
of coefficients in the preceding matrix can be reduced to an equation of 
two unknowns } namely n and nt represented by the primed coeffi- 
cients in the second row of the following primed matrix. This new 
equation and the equation represented by the third row of coefficients 
in the original matrix can then be reduced to an equation of two un- 
knowns, namely and . represented by the primed coefficients 

' o,n nt ^ n 

in the third row of the primed matrix. By continuing this process ^from 
the upper left corner to the lower right corner, the original matrix can 
be reduced to a new, primed matrix of the following form: 



10 


C 2,l ^3,1. 

C 2_, 2 C 3,2 


2,10 

C 3,10 



1,11 

d 2 ,ll 

d 3,ll 

"111 



C 2,12 

C 3,12 


C 2,36 C 3,36 

C 2,37 C 3,37 

C 2,38 

The right member of each of the 38 equations is altered, also, as 
the progression through the original matrix occurs • Therefore, the 
right side of each equation is treated as a separate matrix (called the 
right matrix in the program) • 

Equating the bottom row of the left, primed matrix with that of the 
right, primed matrix gives 


C 1 cp = w f 
2,38^38 38 


which can be written as 


38 


W 38 

"i r 

"2,38 


After 930 is thus determined, the flux at point 37 can be calculated 
by moving up one row in the primed matrix and writing 

C 2, 37^37 + C l37 q) 38 = W 37 


£TT- 
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which can be rearranged in the form 


<P 


37 


ui _ n t /n 
57 5,57^38 

C! 

2,37 


This process continues to the upper left corner of the primed matrix, 
and thus the flux at each of the 38 points can be calculated. The pre- 
ceding procedure is a modified matrix reduction technique that is a modi- 
fication of the Crout matrix reduction technique. 


After a value of the flux at each of the 38 radial points in the 
reactor is thus computed for an assumed value of {3^ at each radial 

position in the core (i.e., the three inner regions), the radial distri- 
bution of neutrons born out of a particular energy group, designated as 
3 n calculated by 


ft i = vZ, <p n 
p n,k f,n,urn,k 

It should be noted that at points 11 and 21 two values of P n ^ can be 
calculated, one value based on the of one region and the other value 

based on the 2^ of the adjacent region. For a calculation of the power 

of the reactor by the method outlined herein, both values must be com- 
puted at an interface within the core. Thus, there are 33 values of 
Pn,k for each energy group. It so happens that a running total of Pn,k 
at each radial point can be made as the problem proceeds; therefore, a 
large amount of machine storage is not required to store the values of 
P n ^ For each energy group. During the entire cycle of computing through 

all energy groups, the previously assumed values of 3t,k> or those calcu- 
lated from the previous cycle, are used. 

After the flux distribution and corresponding values of are 

computed for the first, or highest, energy group, the program proceeds 
to the solution of the matrix of the second energy group to determine 
the radial flux distribution and the 33 values of Pn^k that are added 

to the 3 n k of the next higher group. This procedure continues until 
the solution of all energy groups is completed. 
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Then, after the summations of P n ^ k , designated as are com " 

puted for each radial point in the core, a value of the reactor power P 
is computed by the following equation: 


P 



2rtr ]? /Arpt^ k 


= A' 


J 31 

Z r k ArP 
r 0 


t,k 


(18) 


where A' is an arbitrary constant that cancels itself in computing the 
ratio of powers to detennine k eff . In computing the power, a point k 

is assumed to be at the mean radius of a cylinder of thickness h m . Thus, 

figure 1 shows that special equations must be derived for computing the 
power generated between points 0 and 1, which are a distance hj/2 
apart, and for the cylinders of thickness Yijjz on either side of an 
interface within the core. These special equations are listed below: 


Between Tq and r-^: 


r ]^t,k 


T ~2 p t,i 



(19) 


(Note; p.f. ]_ is assumed to represent the average P-^ from the center of 
the reactor to point 1.) 


Between 


r i and r-i -i 

1Q± - LJ - 

2 


r k ArP 


t,k 





( 20 ) 


where it is assumed by numerical integration that 



h„ 


Pt^ = 24 (" P t,9 + 5 Pt,10 + 8 P t,ll) 


( 21 ) 
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(Note: p-j. j corresponds to the value calculated from a value of 

for region I. A similar regional subscript notation is used on P t in 

the following equations . ) 

Substituting equation (2l) into equation (20) gives 


13h, 


V^t^k = “ST (_P t,9 + 5 P t ,10 + 8 ( 3 t,I,ll) 


( 22 ) 


Between r-> i and r i : 

hi 


r k^ p t,k = 



(m + ~v) |3 t dr =/( ( 10h i + -r)^ 


(23) 


where 


,r 1 

ll| 

V = ( ® P t, 11,11 + 5P t,12 " P t,13 } 

r ll 

Then, substitution of equation (24) into equation (23) gives 

r t^rP t ,k = 9 ^ ( 40h I + h Il)( 8 p t,II,ll + 5 Pt,12 " Pt,13) 

At points 21 and 31 equations representing r-^ArP^^jj. 
a manner similar to that used at point 11 . 


(24) 


(25) 


can be derived in 
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Then, substituting into equation (l8) the right members of equations 
(19), (22), (25), and the equations for points 21 and 31 and the interior 
points of each region results in the following equation: 

|r - + *1^1 + *&>)) + 




k + *3% + *4*4 + 


1111 2-2 r k P t,III,k + y 5 Z 5 + y 6 Z 6 

) 


where 


Y i = - 


z i ~ p t,i 


Y?> = — 


Z 2 = ~Pt,9 + 5(3 t,10 + 8 Pt,I,ll 


40h + h 


Z 3 ~ 8 &t,II,ll + 58 t,12 “ P t,13 


40h I + 39h 


Z 4 ~ _|3 t,19 + 5P t,20 + 8 ^t,II,21 


Y 5 " 


40hj + 40hjj + hjjj 
96 


Z 5 - Q Pt,III,21 + 5 ^t,22 - H,zs 


40hj + 40hjj + 39hjjj 


Z 6 8 t,29 + 5 Pt,30 + 8 Pt,HI,31 


With the value of P/A' calculated for a previous set of values of 
0. , and a new value of P/A' for the values of 0 . obtained from 

y iv X , iC 

the flux values in the reduced matrix solution, a value of k is de- 
fined as 


(P/A’), 


^eff ,new 7p/ A t 


'previous 
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This value of k eff is then compared with a previous value of k e ff to 
determine whether or not the change in k e ff , noted herein as A is 

less than an assumed value (e.g., 10“ 4 ) that specifies a satisfactory 
limit on the convergence. If the limit on Ak is not met, the new 

values of 0 t k and P/A f are normalized by dividing each one by k eff 
(eq. (27)). The normalized values of are then used to compute a 

new flux distribution that results in a new p/A 1 . A new k e ^ is de- 
termined by equation (27), and the Ak e ff is inspected for being less 
than, for example, 10“ 4 . This inspection and recalculation of the prob- 
lem continues until the limits on Ak^^ are met. A discussion of the 

values punched out during, and at the end of, the problem are discussed 
subsequently. 


CALCULATION PROCEDURE 

The solution for the effective multiplication factor k eff of a 

reactor by means of the equations and procedure indicated in the previous 
section is begun by assuming an initial radial distribution of neutrons 
born and then dividing the neutrons into energy groups on the basis of 
the division of the fission spectrum. On the basis of this distribution 
of neutrons born, a radial and energywise flux distribution is calculated. 
This new flux distribution results in a new radial distribution of neutrons 
born and new values of reactor power and k e ff • During each subsequent 
iteration the radial distribution of neutrons born resulting from the 
previous cycle of flux calculations is used to compute new values of re- 
actor power and This iteration procedure continues until the value 

°f k eff has converged within the limits specified. Then, one more time 

through the problem causes the program to punch out: (l) the radial dis- 

tribution of neutron production as a summation of neutrons born out of 
each energy group plus those born out of all higher groups, (2) a normalized 
flux distribution, and (3) the radius of each radial point and the values of 
k g f f and Ak^^^ . 

The details of the program to perform the solution described are in- 
cluded in the appendixes, and an outline of the machine program is given 
in the form of a block diagram in figure 2 0 Appendix B lists the computing 
machine components used and gives instructions with regard to use of the 
machine program. Appendix C discusses the use of the subject program for 
less than nine energy groups and four, or less, concentric regions. Ap- 
pendix D presents an outline and discussion of the various sections of 
the program and gives a complete listing of the machine program. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, October 1, 1958 
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APPENDIX A 


A,A' 

a 

B z 

C 

C' 

D 

d 

d' 

E 

G 

h 

k eff 

^ff 

n 

P 

r 

W 


W f 


SYMBOLS 

coefficients relating power to total number of neutrons 
born (see eq. (18)) 

coefficient defined by equation (16) 

geometric buckling in axial direction 

flux coefficient defined by equations (9) and (10) used 
in forming matrix 

flux coefficient in reduced matrix 
diffusion coefficient, cm 

flux coefficient defined by equations (16) and (17) 
flux coefficient in reduced matrix 

distance from axis of reactor to the halfway point between 
adjacent numbered points k and k+1 (fig. l(b)), cm 

flux coefficient defined in equation (9) 

thickness of concentric cylinders within a region of the 
reactor (fig. 1(b)), cm 

effective multiplication factor 

change in effective multiplication factor 

neutron energy group (a, b, c, ••• i, in order of decreas- 
ing energy to thermal group i) 

reactor total power defined by equation (18), neutrons/sec 

radius, cm 

term in matrix equations representing the neutrons enter- 
ing a group as defined by equation (9) 

term representing neutrons entering a group in the re- 
duced, or primed, matrix 


Y 1 > Y 2>'”> Y 6 


coefficients defined by equation (26), cm 
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o 


Z- 1 .Z 0 .....Z a terms representing neutrons born in equation (26), 

1 neutrons/ (cm 3 ) (sec) 

a ratio of neutrons born into an energy group to the total 

number born 


neutrons born into any one energy group per unit volume 
per unit time, neutrons/(cm 3 )(sec) 

total number of neutrons born into all energy groups per 
unit volume per unit time, neutrons/ (cm 3 ) (sec) 


v 


average number of neutrons produced per fission 


a 


sum of neutron cross sections representing neutron loss 
from a given volume element other than by diffusion, cm 

neutron absorption cross section, cm"! 

neutron fission cross section, cm -1 

neutron slowing-dovm cross section specifying the number 
of neutrons slowing out of one energy group into the 
next lower group, cm -1 

neutron flux, neutrons/(cm £j ) (sec) 


Subscripts : 
a,b,c » • »i 
k 

m 

0,1, 2, 3,. ..39 
I, II, III, IV 


nine neutron energy groups 

numbered radial point in reactor (0^ 2^ 3, • • • as 

shown in fig. l) 

numbered region (l, II ; III , and IV as shown in fig* l) 
radial points 
radial regions 
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APPENDIX B 


INFORMATION ON USE OF PROGRAM 

The machine for which the subject program was developed, is basically 
an IBM 650 magnetic drum data-processing machine (ref. 2) with an IBM 653 
high-speed storage unit (ref. 3) and additional features including in- 
dexing accumulators A, B, and C and a floating decimal arithmetic unit 
(ref. 4). The language in which the program of appendix D is written 
is an optimal programming method called SOAP II (Symbolic Optimal Assembly 
Program II ) developed for the IBM 650 computer (ref. 5). The results of 
processing the original program by the method of SOAP II into the basic 
machine language is shown in the appendix D tabulation, where the original 
coded program is shown on the left and the program in basic machine 
language is given in the four columns of numbers on the right. These 
four columns of numbers represent, column by column from left to right, 

(a) Location of command 

(b) Machine operation code 

(c) Address of data used with machine operation 

(d) Address of next instruction 

The plus signs associated with each command, or the data, are not indi- 
cated; however, the sign, either plus or minus, should be punched in the 
cards to be used as the program. The manner in which the machine pro- 
gram is to be punched into cards , of course, depends on the wiring of 
the control panel in the card reader unit, which varies from one comput- 
ing machine to another. Therefore, general card punching instructions 
cannot be given and are not included in this report. 

The order of the punched cards to be read into the machine is 

1. Suitable load - punch subroutine 


2. Program of appendix D 

3. Set of input data as specified in table I 

4. Basic load card to transfer control to the first instruction of 

program, which is 0500 

When these four items are read in, the machine will compute and punch 
out answers until the end of the program is reached, where there is a 
HALT command. 


OTT 
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When the first run is completed, a second run with only the changes 
desired can be performed by reading the changes into the computer and 
transferring control to address 0500. 

The input locations are specified in table I. The location of 
each input value is dependent on both the physical composition of the 
region and/or the region location; therefore, each storage location has 
been specified in detail in table I. 

The output locations are specified by the set of punch constants 
(program instructions 837 to 844) used in conjunction with an output (or 
PUNCH) subroutine. The instructions in the program that cause punching 
of the output are shown in table II in the order of occurrence during 
the program. Only the values of k e ff and Ak eff are punched out dur- 
ing each iteration to indicate the progress of the convergence of & e ff 
to the final value punched out at the end of the program. 

The means used for loading and punching out data can be adapted to 
the requirements of a similar IBM 650 machine by writing loading and 
punching commands in the drum storages reserved for this purpose (namely, 
1871 to 1999) and utilizing high-speed storages 9055 to 9059 during the 
punching operation. This would mean a possible revision of the program 
instructions used for this purpose, which are listed in table II, and 
the punch constants (837 to 844) in appendix D. The method used by the 
authors for reading the input and punching the output is not discussed 
herein, because the subject program is adaptable to any method presently 
in use by the reader. The authors 1 method is obtainable from the authors 
or the Machine Computing Branch at the NASA Lewis Research Center. The 
system used by the authors is designed to suit many other projects and 
cannot be recommended as a specific means of performing the functions re- 
quired by this problem. 

The only addresses not used that could be used for other purposes 
are addresses 0000, 0037, and 0332. It is understood that locations 
1871 to 1999 are also available for a load-punch routine, as previously 
discussed. 

The program is written in such a way that 1.0 is the original guess 
for each of the 33 values of (3 in the inner three regions . In some 
cases that have been run, this 0 distribution leads to negative fluxes 
in the first iteration. If this difficulty should arise, steps 373 to 377 
can be bypassed by loading the following correction to replace step 372: 

0330 21 0200 0390 

Then it will be necessary to load with the input data the 33 original 
values of 0 stored in addresses 9019 to 9051. 
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APPENDIX C 


USE OF PROGRAM FOR LESS THAN NINE ENERGY GROUPS 

AND FOUR OR LESS REGIONS " t? 

The subject program can be used for any number of neutron energy 
groups up to nine and any number of regions up to four. If fewer than 
nine energy groups are used, input values for the unused energy groups 
must still be read into the machine as discussed in the next paragraph. 

The reason for this is that the machine program generates input for nine 
matrices regardless of the number of groups specified. However, after 
program step 155 (see appendix D) , the program computes only the number 
of groups desired. 

In a solution of fewer than nine groups, the valid cross-section 
values should be located toward the end of each storage region in which 
the input data are stored (see table I) . For those groups not used, the 
cross-section values of 2 .2 and vZ should be zero. Also, the values 

1 f 

of a for the energy groups not used should be zero. The value of 2^ 

for the first group for each region should also be zero, and these zeros 
have to be read in with the other input; whereas, In the case of nine 
groups, steps 832 to 835 provide the necessary zeros for the four re- 
gions. The values of D for the energy groups not used should be equal 
to 1 because the D values appear in the denominator of expressions 
used In the program, and zeros in this case would cause the computing 
machine to stop. 

With this storage arrangement, in the case of nine energy groups 
it is necessary to change the basic machine language of program steps 
158 and 519, which control the input to the matrix calculations. An 
example for six groups would be as follows: 

Step 158 = 0407 83 0005 0463 

Step 519 = 0257 83 0005 0363 

The underlined numbers are equal to 1 less than the number of groups 
being used. 

If fewer than four regions are desired, it is necessary to dupli- 
cate the cross-section values and diffusion coefficients in adjacent 
regions and reduce the value of h in the adjacent regions. In other 
words, one region of the reactor can be treated as two identical regions 
insofar as the machine program is concerned. Thus, a three-region re- 
actor can be treated as having four regions, and the machine program can 
proceed with its solution lor four regions. 


CTT 
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APPENDIX D 


COMPUTING MACHINE PROGRAM INSTRUCTIONS 

The subject program was written for the computing machine described 
in appendix B. In order that the organization of the machine program 
can be better understood, an outline of the program is given in figure 
2. The numbers in the parentheses in each block of figure 2 specify the 
machine instructions that perform the operations indicated. In the tab- 
ulation of the instructions included in this appendix, these instruction 
numbers appear at the extreme left, (in the discussion in this appendix, 
the numbers in parentheses refer to these same instruction numbers.) 

To utilize the indexing features of the computing machine, data are 
stored in groups of consecutive storage locations known as regions (3 to 
18) o The machine loading and punching subroutines are contained in a 
separate group of instructions that are read into the machine ahead of 
the subject program and stored in locations 1871 to 1999 (reserved by 
instruction 22) . 

As indicated in figure 2, there are parts of the program (3 to 379) 
that are performed at the beginning of the problem and are not repeated 
during each iteration to stabilize The remainder of the program 

(383 to 81l) is repeated during each iteration, as shown by the arrow 
indicating the return of the program to instruction (385). Inside of 
this outer loop is an inner loop in the program that computes the neu- 
trons born out of each energy group (518 to 811). As soon as k e ff has 
been stabilized, the program is set for the last iteration, and punching 
takes place during this iteration (746 to 749, 764, 777, 790, 796 to 805, 
494 to 498) . After the last iteration, the program is stopped, and re- 
starting the problem to compute the next case can be done as outlined in 
appendix B. 

Instruction numbers that are not indicated on the outline in figure 
2 are instructions that are not machine commands, such as heading and 
spacing instructions or constants used in the program. The constants 
used are as follows : 

(815 to 823) Numbers 1 to 9 used for code designation of each 
energy group and as constants 

(825 to 830) Other constants that appear in the equations 

(832 to 835) Zero values for the slowing-in cross section of the 
first energy group in each region of the reactor 
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(837 to 844) Code numbers used with punch subroutine to prescribe 
a certain format for the output 


(846 to 847) Values used in the left matrix 

(849) Limit of convergence of k e ff (in this case, 10“^) 


1 1 NINE GROUP REACTOR PROBLEM 

2 1 


3 

REG A0753 

0761 

ALPHAS 

4 

REG B 1767 

1799 

SUM BETAS 

5 

REG C 1720 

1801 

BDRY D 

6 

REG D0762 

0797 

D COEFFS 

7 

REG E0798 

0835 

E VALUES 

8 

REG F 1 729 

1766 

FLUX 

9 

REG G 1840 

1866 

F OF ALPHA 

10 

REG H0717 

0752 

F OF SIG 0 

11 

REG J0681 

0716 

G VALUES 

12 

REG M0836 

1573 

MATRICES 

13 

REG N 1602 

1839 

NORM FLUX 

14 

REG Q 1574 

1609 

SIGMA G 

15 

REG R 1610 

1647 

PT RADII 

16 

REG S 1648 

1683 

SIGMA 

17 

REG V 1684 

1710 

NU SIGMA F 

18 

REG W1711 

1719 

H SO A 

19 

SYN BEGIN 

0500 

START HERE 

20 

EOU PUNCH 

1930 

RKT PUNCH 

21 

BLR 0000 

0000 

BLR 0000 

22 

BLR 1871 

1999 

LOAO PUNCH 

23 

EOU KCHEK 

9018 

TEST FOR 

24 1 



LAST 

25 1 



ITERATION 

26 1 




27 1 

FUNCTIONS 

OF INPUT 

VALUES 


28 1 


29 BEGIN 

RSA 

0003 


GO 

1 



0500 

81 

0003 

0556 

30 GO 

1 

RAU 

9217 






0556 

60 

9217 

0471 

31 


FMP 

8003 




H 

SQUARE 

0471 

39 

8003 

0325 

32 


STU 

9213 






0325 

21 

9213 

0267 

33 


N2A 



GO 

2 



0267 

40 

0320 

0521 

34 


RAU 

9217 






0320 

60 

9217 

0571 

35 


FDV 

9218 




H 

RATIO 

0571 

34 

9218 

0118 

36 


STU 

9260 






0118 

21 

9260 

0464 

37 


AXA 

0001 


GO 

1 



0464 

50 

0001 

0556 

38 1 

39 GO 

2 

RSB 

0008 




H 

SQ A 

0521 

83 

0008 

0177 

40 


RSC 

0035 


GO 

4 



0177 

89 

0035 

0434 

41 GO 

4 

RAU 

S0036 

C 





0434 

60 

7683 

0437 

42 


FDV 

00036 

c 





0437 

34 

6797 

0447 


E-113 



E-113 


23 


43 




FMP 

9210 





0447 

39 

9210 

0210 

44 




FMP 

MIN 1 





0210 

39 

0466 

0317 

45 




STU 

J 0036 

C 




0317 

21 

6716 

0519 

46 




NZB 



GO 

5 


0519 

42 

0172 

0673 

47 




AXB 

0001 


GO 

3 


0172 

52 

0001 

0178 

48 


GO 

5 

AXA 

0001 





0673 

50 

0001 

0429 

49 




SXB 

0008 





0429 

53 

0008 

0285 

50 




NZC 

GO 3 


GO 

6 


0285 

48 

0178 

0439 

51 


GO 

3 

AXC 

0001 


GO 

4 


0178 

58 

0001 

0434 

52 

1 













53 


GO 

6 

RSB 

0026 


GO 

7 


0439 

83 

0026 

0446 

54 


GO 

7 

RAU 

D 0036 

B 




0446 

60 

4797 

0552 

55 




FOV 

D0027 

B 



D RATIO 

0552 

34 

4788 

0038 

56 




STU 

9456 





0038 

21 

9456 

0359 

57 




NZB 



GO 

8 


0359 

42 

0412 

0514 

58 




AXB 

0001 


GO 

7 


0412 

52 

0001 

0446 

59 

1 













60 


GO 

3 

R SC 

0000 




BDR V COEFS 

0514 

89 

0000 

0370 

61 




RSA 

0002 




FOR LEFT 

0370 

81 

0002 

0226 

62 




SXB 

0008 


GO 

10 

MATRIX 

0226 

53 

0008 

0132 

63 


GO 

10 

RAU 

9630 





0132 

60 

9630 

0335 

64 




FMP 

9259 





0335 

39 

9259 

0409 

65 




STU 

C0001 

C 



DA 08 D12 

0409 

21 

7720 

0174 

66 




NZB 



GO 

11 


0174 

42 

0227 

0228 

67 




AXB 

0001 


GO 

9 


0227 

52 

0001 

0484 

68 


GO 

11 

NZ A 



GO 

12 


0228 

40 

0331 

0182 

69 




AXA 

0001 





0331 

50 

0001 

0487 

70 




SXB 

0008 


GO 

9 


0487 

53 

0008 

0484 

71 


GO 

9 

AXC 

0001 


GO 

10 


0484 

58 

0001 

0132 

72 

1 













73 


GO 

12 

RSC 

0026 


GO 

13 


0182 

89 

0026 

0088 

74 


GO 

13 

RAU 

C 002 7 

C 




0088 

60 

7746 

0602 

75 




FMP 

C0082 




D3 07 Dll 

0602 

39 

1801 

0652 

76 




STU 

C 00 54 

C 




0652 

21 

7773 

0276 

77 




RAU 

C 0027 

C 




0276 

60 

7746 

0304 

78 




FAD 

9001 





0304 

32 

9001 

0534 

79 




FMP 

9093 




D2 D6 DIO 

0534 

39 

9003 

0537 

80 




STU 

C 00 8 1 

C 




0537 

21 

7800 

0354 

81 




NZC 



GO 

14 


0354 

48 

0557 

0158 

82 




AXC 

0001 


GO 

13 


0557 

58 

0001 

0088 

83 

1 













84 


GO 

14 

RSA 

0003 




H SO SIG Q 

0158 

81 

0003 

0564 

85 




RSB 

0008 


GO 

16 

OVER D 

0564 

83 

0008 

0420 

86 


GO 

16 

RAU 

9213 





0420 

60 

9213 

0367 

87 




FMP 

00001 

c 




0367 

39 

7574 

0224 

88 




FDV 

D 000 1 

c 




0224 

34 

6762 

0462 

89 




FMP 

MIN 1 





0462 

39 

0466 

0417 

90 




STU 

H 000 1 

c 




0417 

21 

6717 

0470 

91 




NZB 



GO 

17 


0470 

42 

0274 

0324 

92 




AX9 

0001 


GO 

15 


0274 

52 

0001 

0180 

93 


GO 

17 

NZ A 



GO 

18 


0324 

40 

0277 

0278 

94 




AXA 

0001 





0277 

50 

0001 

0584 

95 




SXB 

0008 


GO 

15 


0584 

53 

0008 

0180 



24 


96 


GO 

15 

AXC 

0001 


GO 

16 



0180 

58 

0001 

0420 

97 

1 














98 


GO 

18 

RSA 

0002 




H SQ 

ALPHA 

0278 

81 

0002 

0634 

99 




RSB 

0008 




OVER 

D 

0634 

83 

0008 

0090 

100 




RSC 

0026 


GO 

20 



0090 

89 

0026 

0496 

101 


GO 

20 

RAU 

A0009 

B 





0496 

60 

4761 

0467 

102 




N2U 



GO 

21 



0467 

44 

0621 

0222 

103 




FDV 

00027 

C 





0621 

34 

6788 

0138 

104 




FMP 

MIN 1 






0138 

39 

0466 

0517 

105 




FMP 

9212 


GO 

21 



0517 

39 

9212 

0222 

106 


GO 

21 

STU 

G0027 

C 





0222 

21 

7866 

0569 

107 




NZB 



GO 

22 



0569 

42 

0272 

0374 

108 




AXB 

0001 


GO 

19 



0272 

52 

0001 

0328 

109 


GO 

22 

NZA 



GO 

23 



0374 

40 

0327 

0378 

110 




AX A 

0001 






0327 

50 

0001 

0385 

711 




SXB 

0008 


GO 

19 



0385 

53 

0008 

0328 

112 


GO 

19 

AXC 

0001 


GO 

20 



0328 

58 

0001 

0496 

113 

1 














114 


GO 

23 

RSA 

0003 




R AND 

E 

0378 

81 

0003 

0435 

115 




RSB 

0008 




CALC 


0435 

83 

0008 

0641 

116 




RSC 

0028 






0641 

89 

0028 

0497 

117 




RAU 

9014 






0497 

60 

9014 

0505 

118 




FDV 

9002 






0505 

34 

9002 

0208 

119 




STU 

R 000 1 






0208 

21 

1610 

0614 

120 




STU 

9000 






0614 

21 

9000 

0671 

121 




STU 

E0001 






0671 

21 

0798 

0404 

122 




FAD 

9000 






0404 

32 

9000 

0485 

123 




STU 

R0002 


GO 

24 



0485 

21 

1611 

0664 

124 


GO 

24 

FAD 

9000 






0664 

32 

9000 

0444 

125 




STU 

E0030 

c 





0444 

21 

6827 

0230 

126 




FAD 

9000 






0230 

32 

9000 

0459 

127 




STU 

R 003 1 

c 





0459 

21 

7640 

0494 

128 




NZB 



GO 

25 



0494 

42 

0547 

0498 

129 




AXB 

0001 


GO 

26 



0547 

52 

0001 

0454 

130 


GO 

26 

AXC 

0001 


GO 

24 



0454 

58 

0001 

0664 

131 


GO 

25 

NZA 



JCALC 



0498 

40 

0504 

0554 

132 




AXA 

0001 






0504 

50 

0001 

0260 

133 




RAU 

9217 






0260 

60 

9217 

0322 

134 




FDV 

9002 






0322 

34 

9002 

0375 

135 




STU 

9000 






0375 

21 

9000 

0535 

136 




RAU 

R0031 

c 





0535 

60 

7640 

0546 

137 




NZC 



GO 

27 



0546 

46 

0349 

0604 

138 




SXB 

0009 


GO 

26 



0349 

53 

0009 

0454 

139 


GO 

27 

SXB 

0006 


GO 

26 



0604 

53 

0006 

0454 

140 

1 














141 


JCALC 

RSA 

0008 


GO 

90 

J VALUES 

0554 

81 

0008 

0310 

142 


GO 

90 

RAU 

J0009 

A 





0310 

60 

2689 

0544 

143 




STU 

W0009 

A 





0544 

21 

3719 

0372 

144 




FSB 

9002 






0372 

33 

9002 

0654 

145 




STU 

J0009 

A 





0654 

21 

2689 

0292 

1 46 




NZA 



GO 

91 



0292 

40 

0596 

0646 

147 




AXA 

0001 


GO 

90 



0596 

50 

0001 

0310 

148 


GO 

91 

RSA 

0026 


GO 

92 



0646 

81 

0026 

0555 


E-113 
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149 

150 

151 

152 

153 

154 1 

155 1 

156 1 

157 

158 
1 5 9 
160 
161 
162 
16 3 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 
1 76 

177 

178 

179 1 

180 1 
181 1 


GO 92 


RAU J0036 A 

FSB 9002 

STU J0036 A 

N2 A CROUT 

AX A 0001 GO 92 

LEFT MATRIX POINTS 1-10 


CROUT 


GO 34 


MATRIX 
GROUP 
NET POINT 


GO 28 


RAC 0000 

RS8 0008 GO 34 

RSA 0008 

RAU W0009 B 

FDV 9002 

FSB 9003 

STU M0001 C 

RAU 9003 

FDV M0001 C 

STU M0041 C GO 28 

FMP E0009 A 

STU 9000 

RAU R0010 A 

FMP J 0009 B 

FSB 9000 

STU M0002 C 

RAU E0010 A 

FDV MQ002 C 

STU M 0042 C 

NZA BDRY1 

AXA 0001 

AXC 0001 GO 28 


LEFT MATRIX POINTS 11 AND 12 


182 BDRY1 

sxc 

0008 

183 

RAU 

C0082 

184 

FSB 

M0049 

185 

STU 

M0081 

186 

FMP 

M0050 

187 

STU 

9000 

188 

RAU 

C0063 

189 

FSB 

9000 

190 

STU 

MOO 1 1 

191 

RAU 

C0036 

192 

FDV 

MOO 1 1 

193 

STU 

M0051 

194 

RAU 

C 0009 

195 

FDV 

M0011 

196 

STU 

M 0 0 8 2 

197 1 



198 

RAU 

E0011 

199 

FMP 

MOOS 1 

200 

STU 

9000 

201 

RAU 

R 00 1 2 


POINT 11 


C 

C 

c 


c 

B 

C 

C 

B 

C 

C 


POINT 12 


0555 60 

0422 33 

0605 21 

0619 40 

0472 50 


0251 88 

0257 83 

0363 81 

0369 60 

0373 34 

0176 33 

0305 21 

0189 60 

0297 34 

0136 21 

0279 39 

0356 21 

0413 60 

0423 39 

0239 33 

0419 21 

0040 60 

0311 34 

0287 21 

0130 40 

0383 50 

0289 58 


0284 59 

0140 60 

0655 33 

0511 21 

0669 39 

0585 21 

0594 60 

0587 33 

0567 21 

0399 60 

0509 34 

0597 21 

0489 60 

0635 34 

0647 21 

0520 60 

0617 39 

0336 21 

0644 60 


2716 0422 
9002 0605 
2716 0619 
0472 0251 
0001 0555 


0000 0257 
0008 0363 
0008 0369 
5719 0373 

9002 0176 

9003 0305 
6836 0189 
9003 0297 

6836 0136 

6876 0279 

2806 0356 
9000 0413 
3619 0423 
4689 0239 
9000 0419 

6837 0040 

2807 0311 
6837 0287 

6877 0130 
0383 0284 

0001 0289 
0001 0279 


0008 0140 
1801 0655 

6884 0511 

6916 0669 

6885 0585 
9000 0594 
5782 0587 
9000 0567 
6846 0399 
5755 0509 
6846 0597 

6886 0489 
5728 0635 
6846 0647 

6917 0520 

0808 0617 
6886 0336 
9000 0644 
1621 0425 
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2 0 2 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 : 

214 ; 

215 ] 

216 
217 
2 18 

219 

220 
221 
222 

223 

224 

225 

226 
227 
2 28 

229 1 

230 1 

231 1 

232 

233 

234 

235 

236 
2 37 

238 

239 

240 

241 

242 

243 
2 44 

245 

246 

247 1 

248 

249 

250 

251 

252 

253 

254 


GO 29 


BDRY2 


FMP 

J0018 B 



0425 

39 

4698 

0548 

FSB 

9000 



0548 

33 

9000 

0377 

STU 

MOO 1 2 C 



0377 

21 

6847 

0606 

RAU 

E 00 1 1 



0606 

60 

0008 

0667 

FMP 

MOOS 2 C 



0667 

39 

6917 

1867 

STU 

9000 



1867 

21 

9000 

0475 

RAU 

E0012 



0475 

60 

0809 

0168 

FSB 

9000 



0168 

33 

9000 

0598 

FDV 

M0012 C 



0598 

34 

6847 

0648 

STU 

M0052 C 



0648 

21 

6887 

0190 

SXA 

0007 

GO 

29 

0190 

51 

0007 

0301 

LEFT MATRIX 

POINTS 13-20 





RAU 

E0019 A 



0301 

60 

2816 

0221 

FMP 

MOOS 2 C 



0221 

39 

6887 

0337 

STU 

9000 



0337 

21 

9000 

0595 

RAU 

R 0020 A 



0595 

60 

3629 

0433 

FMP 

J0018 B 



0433 

39 

4698 

0348 

FSB 

9000 



0348 

33 

9000 

0127 

STU 

MOO 1 3 C 



0127 

21 

6848 

0351 

RAU 

E 002 0 A 



0351 

60 

2817 

0271 

FDV 

MOO 1 3 C 



0271 

34 

6848 

0398 

STU 

MOOS 3 C 



0398 

21 

6888 

0541 

NZA 


BDRY2 

0541 

40 

0394 

0645 

AXA 

0001 



0394 

50 

0001 

0401 

AXC 

0001 

GO 

29 

0401 

58 

0001 

0301 


LEFT MATRIX POINTS 21 AND 22 


sxc 

0007 


POINT 21 

0645 

59 

0007 

0656 

RAU 

C0082 



0656 

60 

1801 

0607 

FSB 

M0059 

C 


0607 

33 

6894 

0522 

STU 

M0079 

C 


0522 

21 

6914 

0218 

FMP 

M0060 

C 


0216 

39 

6895 

0449 

STU 

9000 



0449 

21 

9000 

0657 

RAU 

C0072 

B 


0657 

60 

5791 

0499 

FSB 

9000 



0499 

33 

9000 

0479 

STU 

M0021 

C 


0479 

21 

6856 

0559 

RAU 

C 0045 

B 


0559 

60 

5764 

1869 

FDV 

MOO 2 1 

C 


1869 

34 

6856 

0258 

STU 

M 006 1 

C 


0258 

21 

6896 

0549 

RAU 

COO 1 8 

B 


0549 

60 

5737 

0342 

FDV 

M 002 1 

C 


0342 

34 

6856 

0308 

STU 

M 0039 

c 


0308 

21 

6874 

0427 


RAU E0021 
FMP M0061 C 
STU 9000 
RAU R0022 
FMP J0027 B 
FSB 9000 
STU M0022 C 


POINT 22 


0427 60 0818 0424 
0424 39 6896 0599 
0599 21 9000 0358 
0358 60 1631 0386 
0386 39 4707 0408 
0408 33 9000 0637 
0637 21 6857 0360 



CY-4 back 


27 


H 

I 

A 


2 55 




RAU 

E 002 1 






0360 

60 

0818 

0474 

256 




FMP 

M0039 

C 





0474 

39 

6874 

0524 

257 




STU 

9000 






0524 

21 

9000 

0381 

258 




RAU 

E 0022 






0381 

60 

0619 

0574 

259 




FSB 

9000 






0574 

33 

9000 

0458 

260 




FDV 

M0022 

C 





0458 

34 

6857 

0508 

261 




STU 

M01'62 

C 





0508 

21 

6897 

0558 

262 




SX A 

0007 


GO 

30 



0558 

51 

0007 

0451 

263 

1 














264 

1 



LEFT MATRIX 

POINTS 

23-3 0 






265 

1 














266 


GO 

30 

RAU 

E 0029 

A 





0451 

60 

2826 

0231 

267 




FMP 

M0062 

r 





0231 

39 

6897 

0347 

268 




STU 

9000 






0347 

21 

9000 

0355 

269 




RAU 

R 0030 

A 





0355 

60 

3639 

0393 

270 




FMP 

J0027 

B 





0393 

39 

4707 

0307 

271 




FSB 

9000 






0307 

33 

9000 

0387 

272 




STU 

MOO 2 3 

C 





0387 

21 

6858 

0361 

273 




RAU 

E0030 

A 





0361 

60 

2827 

0281 

274 




FDV 

M0023 

C 





0281 

34 

6858 

0058 

275 




STU 

M 006 3 

c 





0058 

21 

6898 

0501 

276 




NZ A 



GO 

31 



0501 

40 

0054 

0405 

277 




AX A 

0001 






0054 

50 

0001 

0110 

278 




AX C 

0001 


GO 

30 



0110 

58 

0001 

0451 

279 

1 














280 

1 



LEFT MATRIX 

POINTS 

31 AND 

32 





281 

1 














282 


GO 

31 

sxc 

0007 




POINT 

31 

0405 

59 

0007 

0561 

283 




RAU 

C 0 0 8 2 






0561 

60 

1801 

0608 

284 




FSB 

M 0069 

c 





0608 

33 

6904 

0431 

285 




STU 

M0080 

c 





0431 

21 

6915 

0268 

286 




FMP 

M0070 

c 





0268 

39 

6905 

0658 

287 




STU 

9000 






0658 

21 

9000 

0318 

288 




RAU 

CO 081 

B 





0318 

60 

5800 

0609 

289 




FSB 

9000 






0609 

33 

9000 

0539 

290 




STU 

MOO 3 1 

C 





0539 

21 

6866 

0570 

291 




RAU 

C 0054 

B 





0570 

60 

5773 

0477 

292 




FDV 

M0031 

c 





0477 

34 

© 866 

0368 

293 




STU 

M0071 

c 





0368 

21 

6906 

0659 

294 




RAU 

C 0027 

B 





0659 

60 

5746 

0410 

295 




FDV 

M 003 1 

C 





0410 

34 

6866 

0418 

296 




STU 

M0040 

c 





0418 

21 

6875 

0428 

297 

1 














298 




RAU 

E 003 1 




POINT 

32 

0428 

60 

0828 

0436 

299 




FMP 

M0071 

c 





0436 

39 

6906 

0460 

300 




STU 

9000 






0460 

21 

9000 

0468 

301 




RAU 

R 0032 






0468 

60 

1641 

0649 

302 




FMP 

J0036 

B 





0649 

39 

4716 

0518 

303 




FSB 

9000 






0518 

33 

9000 

0510 

304 




STU 

M0032 

c 





0510 

21 

6867 

0620 

305 




RAU 

E0031 






0620 

60 

0828 

0486 

?06 




FMP 

M 0040 

c 





0486 

39 

6875 

0525 

307 




STU 

9000 






0525 

21 

9000 

0536 



28 


308 




RAU 

E0032 




0536 

60 

0829 

0586 

309 




FSB 

9000 




0586 

33 

9000 

0568 

310 




FDV 

M0032 

C 



0568 

34 

6867 

0618 

3H 




STU 

M0072 

C 



0618 

21 

6907 

0560 

312 




SX A 

0005 


GO 32 


0560 

51 

0005 

0551 

313 

1 












314 

1 



LEFT MATRIX POINTS 

33-38 





315 

1 












316 


GO 

32 

RAU 

E 0037 

A 



0551 

60 

2834 

0339 

317 




FMP 

M0072 

C 



0339 

39 

6907 

0357 

318 




STU 

9000 




0357 

21 

9000 

0415 

319 




RAU 

R0038 

A 



0415 

60 

3647 

0601 

320 




FMP 

J0036 

B 



0601 

39 

4716 

0366 

321 




FSB 

9000 




0366 

33 

9000 

0246 

322 




STU 

M0033 

C 



0246 

21 

6868 

0321 

323 




RAU 

E 0038 

A 



0321 

60 

2 835 

0389 

324 




FDV 

M0033 

C 



0389 

34 

6868 

0068 

325 




STU 

M0073 

c 



0068 

21 

6908 

0411 

326 




NZA 



GO 33 


0411 

40 

0214 

0465 

327 




AXA 

0001 




0214 

50 

0001 

0270 

328 




AXC 

0001 


GO 32 


0270 

58 

0001 

0551 

329 

1 












330 

1 



SET-UP 

MATRIX FOR N 

EXT GROUP 





331 

1 












? 32 


GO 

33 

NZB 



GO 35 


0465 

42 

0668 

0670 

333 




AXB 

0001 




0668 

52 

0001 

0624 

334 




AXC 

0077 


GO 34 


0624 

58 

0077 

0363 

335 

1 












336 

1 



STORE ORIGINAL KEFF 

AND POWER 





337 

1 












338 


GO 

35 

RAU 

9001 




0670 

60 

9001 

0527 

339 




STU 

9011 



ORIG KEFF 

0527 

21 

9011 

0288 

340 




STU 

9010 



ORIG POWER 

0288 

21 

9010 

0376 

341 

1 












342 

1 



COEFFS 

FOR COMPUTING POWER 





343 

1 












344 




RAU 

9014 




0376 

60 

9014 

0338 

345 




FDV 

9008 




0338 

34 

9008 

0442 

346 




STU 

VI 




0442 

21 

0028 

0531 

347 




RAU 

C T 39 




0531 

60 

0388 

0426 

348 




FDV 

CT 12 




0426 

34 

0579 

0629 

349 




STU 

y 2 




0629 

21 

0046 

0476 

350 




RAU 

CT 40 




0476 

60 

0679 

0438 

351 




FMP 

9014 




0438 

39 

9014 

0492 

352 




FAD 

9015 




0492 

32 

9015 

0526 

353 




STU 

9000 




0526 

21 

9000 

0488 

354 




FDV 

CT 96 




0488 

34 

0542 

0592 

355 




STU 

Y 3 




0592 

21 

0066 

0576 

356 




RAU 

CT 38 




0576 

60 

0280 

0538 

357 




FMP 

9015 




0538 

39 

9015 

0642 

358 




FAD 

9000 




0642 

32 

9000 

0626 

359 




STU 

9000 




0626 

21 

9000 

0588 

360 




FDV 

CT 96 




0588 

34 

0542 

0676 



29 


361 



STU 

Y 4 



362 



RAU 

9015 



363 



FAD 

9000 



364 



FAD 

9016 



365 



STU 

9000 



366 



FDV 

C T 96 



367 



STU 

Y 5 



368 



RAU 

C T 36 



369 



FMP 

9016 



370 



FAD 

9000 



371 



FDV 

C T 96 



372 



STU 

Y 6 



373 



RS A 

0032 



374 



RAU 

9001 


MBET A 

375 


MBET A 

STD 

9251 



376 



NZ A 



MKCHK 

377 



AX A 

0001 


MBET A 

378 


MKCHK 

RAU 

MIN 1 



379 



STU 

KCHEK 


ZEROB 

380 

1 






381 

1 


MOVE BETAS ZERO 

382 

1 






383 


MOV EB 

SET 

9019 



384 



LDB 

B 000 1 


ZEROB 

385 


ZEROB 

RS A 

0032 



386 



RAU 

H 000 1 


RESET 

387 


RESET 

STU 

B 003 3 

A 


388 



NZ A 



D OW E R 

389 



A X A 

0001 


RESET 

390 

1 






391 

1 


POWr 9 i 

COMPUTATION 

392 

1 






393 


POWER 

SX A 

0008 



394 



RAU 

9029 



3 9 c 



FMP 

9008 



396 



FSB 

9027 



? 97 



STU 

9000 



398 



RAU 

9028 



3 99 

1 






400 




9005 



4 C 1 



FAD 

9000 



4 O 2 



FMP 

Y 2 



u j 3 



FAD 

9 019 



4 0 4 



FMP 

Y 1 



405 



STU 

9000 


GO 38 

4 0 6 


GO 3 6 

RAU 

RQ010 

A 


<10 7 



FMP 

9228 



40 B 



FAD 

9000 



4 0 ° 



NZ A 



GO 3 9 

4 1C 



STU 

9000 



411 



A X A 

0001 


GO 38 

4 1 ? 


GO 39 

FMP 

9014 



413 



STU 

9000 




0676 21 0064 0577 

0577 60 9015 0638 

0638 32 9000 0627 

0627 32 9016 0677 

0677 21 9000 0509 

0589 34 0542 0528 

0528 21 0102 0578 

0578 60 0280 0639 

0639 39 9016 0628 

0628 32 9000 0678 

0678 34 0542 0330 

0330 21 0200 0380 

038C 81 0032 0290 

0 29C 60 9001 0A30 

0430 24 9251 0480 

0480 40 0340 0390 

0340 50 0001 0430 

0390 60 0466 0530 

0530 21 9018 0572 

3ETA SUMS 

0448 27 9019 0610 

0610 09 1767 0572 

0572 81 0032 0478 

0478 60 0717 0622 

0622 21 3799 0660 

0660 40 1868 0150 

1868 50 0001 0622 


REGION I 0150 51 0008 0156 

0156 60 9029 0063 

0063 39 9008 0016 

0016 33 9027 0295 

0295 21 9000 0253 

0253 60 9028 0011 

0011 39 9005 0014 

0014 32 9000 0093 

0093 39 0046 0096 

0096 32 9013 0075 

0075 39 0028 0078 

0078 21 9000 0135 

0135 60 3619 0023 

0023 39 9228 0026 

0026 32 9000 0105 

0105 0008 0109 

0008 21 9000 0115 

0115 50 0001 0135 

0109 39 9014 0012 

001? 21 9000 0069 



30 


4l4 1 


415 


SX A 

0008 


REGION 

II 

0069 

51 

0008 

0125 

416 


RAU 

9030 




0125 

60 

9030 

0133 

417 


FMP 

9008 




0133 

39 

9008 

0086 

418 


FSB 

9032 




0086 

33 

9032 

0165 

419 


STU 

9052 




0165 

21 

9052 

0073 

420 


RAU 

9031 




0073 

60 

9031 

0031 

421 


FMP 

9005 




0031 

39 

9005 

0084 

422 


FAD 

9052 




0084 

32 

9052 

0113 

'•23 


FMP 

Y 3 




0113 

39 

0066 

0116 

424 


STU 

9052 




0116 

21 

9052 

0123 

425 


RAU 

9040 




0123 

60 

9040 

0081 

426 


FMP 

9008 




0081 

39 

9008 

0134 

427 


FSB 

9038 




0134 

33 

9038 

0163 

428 


STU 

9053 




0163 

21 

9053 

0021 

429 


RAU 

9039 




0021 

60 

9039 

0029 

430 


FMP 

9005 




0029 

39 

9005 

0082 

431 


FAD 

9053 




0082 

32 

9053 

0061 

432 


FMP 

Y 4 




0061 

39 

0064 

0114 

433 


FAD 

9052 




0114 

32 

9052 

0143 

434 


STU 

9052 GO 

40 



0143 

21 

9052 

0051 

435 GO 

40 

RAU 

R 0020 A 




0051 

60 

3629 

0183 

436 


FMP 

9239 




0183 

39 

9239 

0087 

437 


FAD 

9052 




0087 

32 

9052 

0017 

4 38 


NZA 

GO 

41 



0017 

40 

0070 

0071 

439 


STU 

9052 




0070 

21 

9052 

0027 

440 


AXA 

0001 GO 

40 



0027 

50 

0001 

0051 

441 GO 

41 

FMP 

9015 




0071 

39 

9015 

0024 

442 


FAD 

9000 




0024 

32 

9000 

0303 

443 


STU 

9000 




0303 

21 

9000 

0111 

444 1 











445 


RAU 

9041 


REGION 

III 

0111 

60 

9041 

0119 

446 


FMP 

9008 




0119 

39 

9008 

0022 

447 


FSB 

9043 




0022 

33 

9043 

0101 

448 


STU 

9052 




0101 

21 

9052 

0159 

449 


RAU 

9042 




0159 

60 

9042 

0067 

450 


FMP 

9005 




0067 

39 

9005 

0120 

451 


FAD 

9052 




0120 

32 

9052 

0099 

452 


FMP 

Y 5 




0099 

39 

0102 

0152 

453 


STU 

9052 




0152 

21 

9052 

0209 

454 


RAU 

9051 




0209 

60 

9051 

0117 

455 


FMP 

9008 




0117 

39 

9008 

0170 

4 56 


FSB 

9049 




0170 

33 

9049 

0149 

457 


STU 

9053 




0149 

21 

9053 

0007 

458 


RAU 

9050 




0007 

60 

9050 

0215 

459 


FMP 

9005 




0215 

39 

9005 

0018 

460 


FAD 

9053 




0018 

32 

9053 

0097 

461 


FMP 

Y 6 




0097 

39 

0200 

0250 

'•62 


FAD 

9052 




0250 

32 

9052 

0079 

46 3 


STU 

9052 




0079 

21 

9052 

0137 

464 


SX A 

0008 GO 

42 



0137 

51 

0008 

0193 

465 GO 

42 

RAU 

R 003 0 A 




0193 

60 

3639 

0243 

466 


FMP 

9250 




0243 

39 

9250 

0298 


sit - a 



E-113 


31 


467 



FAD 

9052 



0298 

32 

9052 

0077 

468 



NZ A 


GO 43 


0077 

40 

0080 

0131 

46 9 



STU 

9052 



0080 

21 

9052 

0187 

470 



AXA 

0001 

GO 42 


0187 

50 

0001 

0193 

471 


GO 43 

FMP 

9016 



0131 

39 

9016 

0636 

472 



FAD 

9000 


NEW POWER 

0636 

32 

9000 

0282 

473 



STU 

9000 



0282 

21 

9000 

0672 

4 7 4 

1 










4 7 5 

1 


NEW KEFF 

AND DEL KEFF PLUS TST 





476 

1 










477 



FDV 

9010 



0672 

34 

9010 

0575 

478 



STU 

9012 


NEW KEFF 

0575 

21 

9012 

0188 

479 



FSB 

9011 



0188 

33 

9011 

0674 

460 



STU 

9013 


DELTA KEFF 

0674 

21 

9013 

0481 

481 



RAU 

KCHEK 



0481 

60 

9018 

0392 

482 



BMI 


R PCH 


0392 

46 

0611 

0661 

433 



RAU 

DEL K 



0611 

60 

0625 

0529 

484 



FSM 

9013 



0529 

38 

9013 

0512 

4 8 5 



STU 

KCHEK 

K PCH 


0512 

21 

9018 

0240 

4 6 6 

1 










467 

1 


PUNCH KEFF AND DEL 

KEFF 





*■8 8 

1 










489 


K PCH 

RAL 

PCKEF 



0240 

65 

0662 

0326 

490 



LDD 

MKEFF 

PUNCH 


0326 

69 

0562 

1930 

491 

1 




\ 






492 

1 


PUNCH R 

KEFF AND 

DEL KEFF 





4 9 3 

1 










4 9 4 


R PCH 

RAL 

PCR 


PCH RAO 1 1 

0661 

65 

0581 

0590 

4 9 5 



LDD 


PUNCH 


0590 

69 

0631 

1930 

4 9 6 



RAL 

PCKEF 


PUNCH KEFF 

0631 

65 

0662 

0232 

4 0 7 



LDD 


PUNCH 

AND DEL K 

0232 

69 

0640 

1930 

4 9 3 
4 '* Q 



HLT 

9999 



0640 

01 

9999 

9999 

5 •. * 0 

U "j 

1 


end of calculation 






c'i , 2 

l 

1 


MOVE KEFF 

AND NORMALIZE POWER 





5 0 4 

l 

MKEFF 

RAU 

9012 



0562 

60 

9012 

0675 

5 ^ 5 



STU 

9011 



0675 

21 

9011 

0238 

5 06 



R S A 

0032 

NORMB 


0238 

81 

0032 

0382 

5 0 7 


NO \'-i3 

RAU 

9251 



0382 

60 

9251 

0432 

3 .j j 



FDV 

9012 



0432 

34 

9012 

0482 

509 



STU 

9251 



0482 

21 

9251 

0532 

510 



NZA 


NORMP 


0532 

40 

0582 

0632 

511 



AXA 

0001 

NORMB 


0582 

50 

0001 

0382 

5 12 


NORM? 

RAU 

9000 



0632 

60 

9000 

0612 

513 



FDV 

9012 



0612 

34 

9012 

1870 

514 



STU 

9010 

COMPW 


1870 

21 

9010 

0651 

515 

1 










516 

1 


RIGHT MATRIX POINT 

1 





517 

1 










518 


COMP W 

RSC 

0000 


MATRIX 

0651 

89 

0000 

0407 

519 



RSB 

0008 

GO 78 

GROUP NO 

0407 

83 

0008 

0463 



32 


520 


GO 78 

RSA 

0008 


NET POINT 

0463 

81 

0008 

0469 

5 21 



RAU 

G0009 

B 


0469 

60 

5848 

0553 

522 



NZU 


GO 

44 

0553 

44 

0457 

0108 

523 



FDV 

9002 



0457 

34 

9002 

0160 

524 



FMP 

9019 

GO 

44 

0160 

39 

9019 

0108 

525 


GO 44 

STU 

9000 



0108 

21 

9000 

0515 

526 



RAU 

H 0009 

B 


0515 

60 

4725 

0329 

527 



NZU 


GO 

45 

0329 

44 

0483 

0334 

528 



FDV 

9002 



0483 

34 

9002 

0186 

5 2 9 



FMP 

F 000 1 

GO 

45 

0186 

39 

1729 

0334 

530 


GO 45 

FAD 

9000 



0334 

32 

9000 

0513 

531 



FDV 

M0001 

C 


0513 

34 

6336 

0236 

532 



STU 

F 0001 

GO 

48 

0236 

21 

1729 

0206 

533 


GO 48 

FMP 

E0009 

A 


0206 

39 

2806 

0406 

534 



STU 

9054 



0406 

21 

9054 

0563 

535 



RAU 

G0009 

B 


0563 

60 

5848 

0603 

536 



NZU 


GO 

46 

0603 

44 

0507 

0300 

537 



FMP 

9228 

GO 

46 

0507 

39 

9228 

0300 

538 

1 










539 

1 


RIGHT MATRIX 

POINTS 2-10 





540 

1 










541 


GO 46 

STU 

9000 



0300 

21 

9000 

0057 

542 



RAU 

H 0009 

B 


0057 

60 

4725 

0129 

543 



NZU 


GO 

47 

0129 

44 

0233 

0184 

544 



FMP 

F0010 

A GO 

47 

0233 

39 

3738 

0184 

545 


GO 47 

FAD 

9000 



0184 

32 

9000 

0213 

546 



FMP 

R0010 

A 


0213 

39 

3619 

0169 

547 



FSB 

9054 



0169 

33 

9054 

0199 

548 



FDV 

M 0002 

C 


0199 

34 

6837 

0237 

549 



STU 

F0010 

A 


0237 

21 

3738 

0391 

550 



N Z A 


GO 

49 

0391 

40 

0294 

0345 

551 



AXC 

0001 



0294 

58 

0001 

0350 

552 



AX A 

0001 

GO 

48 

0350 

50 

0001 

0206 

c 5 3 

1 










5 54 

1 


RIGHT MATRIX 

; POINT 11 





555 

1 










556 


GO 49 

sxc 

0008 



0345 

59 

0008 

0302 

557 



FMP 

M0081 

C 


0302 

39 

6916 

0416 

558 



FAD 

F 0009 



0416 

32 

1737 

0613 

559 



FMP 

MIN 1 



0613 

39 

0 466 

0516 

560 



FDV 

M0011 

C 


0516 

34 

6846 

0296 

561 



STU 

F0011 



0296 

21 

1739 

0142 

562 



SX A 

0008 

GO 

52 

0142 

51 

0008 

0400 

563 

1 










564 

1 


RIGHT 

MATRIX POINTS 12-20 





565 

1 










566 


GO 52 

FMP 

E0019 

A 


0400 

39 

2816 

0166 

567 



STU 

905A 



0166 

21 

9054 

0173 

568 



RAU 

GOO 1 8 

B 


0173 

60 

5857 

0161 

569 



NZU 


GO 

50 

0161 

44 

0265 

0216 

570 



FMP 

9239 

GO 

50 

0265 

39 

9239 

0216 

571 


GO 50 

STU 

9000 



0216 

21 

9000 

0223 

572 



RAU 

H0018 

B 


0223 

60 

4734 

0089 


E- 113 



CY-5 E-113 


33 


573 




NZU 



GO 

51 

0089 

44 

0293 

0344 

574 




FMP 

F 002 0 

A 

GO 

51 

0293 

39 

3748 

0344 

575 


GO 

51 

FAD 

9000 




0344 

32 

9000 

0273 

576 




FMP 

R 0020 

A 



0273 

39 

3629 

0179 

577 




FSB 

9054 




0179 

33 

9054 

0259 

578 




FDV 

MOO 1 2 

C 



0259 

34 

6847 

0147 

579 




STU 

F 0020 

A 



0147 

21 

3748 

0151 

580 




NZA 



GO 

53 

0151 

40 

0004 

0155 

581 




AXC 

0001 




0004 

58 

0001 

0060 

582 




AX A 

0001 


GO 

52 

0060 

50 

0001 

0400 

583 

1 












584 

1 



RIGHT MATRIX 

POINT 21 





585 

1 












586 


GO 

53 

sxc 

0008 




0155 

59 

0008 

0461 

587 




FMP 

M0079 

C 



0461 

39 

6914 

0264 

588 




FAD 

F0019 




0264 

32 

1747 

0473 

589 




FMP 

MIN 1 




0473 

39 

0466 

0566 

590 




FDV 

M0021 

C 



0566 

34 

6856 

0456 

591 




STU 

F 002 1 




0456 

21 

1749 

0352 

592 




SX A 

0008 


GO 

54 

0352 

51 

0008 

0450 

593 

1 












594 

1 



RIGHT 1 

MATRIX 

POINTS 22-30 





595 

1 












596 


GO 

54 

FMP 

E 0029 

A 



0450 

39 

2826 

0076 

597 




STU 

9054 




0076 

21 

9054 

0283 

598 




RAU 

G 002 7 

B 



0283 

60 

5866 

0121 

599 




NZU 



GO 

55 

0121 

44 

0175 

0126 

600 




FMP 

9250 


GO 

55 

0175 

39 

9250 

0126 

601 


GO 

55 

STU 

9000 




0126 

21 

9000 

0333 

602 




RAU 

H0027 

e 



0333 

60 

4743 

0197 

603 




NZU 



GO 

56 

0197 

44 

0201 

0202 

604 




FMP 

F 00 30 

A 

GO 

56 

0201 

39 

3758 

0202 

605 


GO 

56 

FAD 

9000 




0202 

32 

9000 

0181 

606 




FMP 

R 0030 

A 



0181 

39 

3639 

0139 

607 




FSB 

9054 




0139 

33 

9054 

0219 

608 




FDV 

M 002 2 

C 



0219 

34 

6857 

0107 

609 




STU 

F 003 0 

A 



0107 

21 

3758 

0211 

6 1 G 




NZA 



GO 

57 

0211 

40 

0164 

0315 

6 1 1 




AXC 

0001 




0164 

58 

0001 

0220 

612 




AXA 

0001 


GO 

54 

0220 

50 

0001 

0450 

613 

1 












6 1 4 

1 



RIGHT 

MATRIX 

POINT 31 





615 

1 












6 1 6 


GO 

57 

sxc 

0008 




0315 

59 

0008 

0371 

6 1 7 




FMP 

M 008 0 

C 



0371 

39 

6915 

0565 

6 18 




FAD 

F 0029 




0565 

32 

1757 

0533 

6 19 




FMP 

MIN 1 




0533 

39 

0466 

0616 

620 




FDV 

M0031 

C 



0616 

34 

6866 

0666 

621 




STU 

F 003 1 




0666 

21 

1759 

0262 

622 




SX A 

0006 


GO 

58 

0262 

51 

0006 

0650 


623 1 

624 1 

625 1 


RIGHT MATRIX POINTS 32~38 



34 


626 GO 58 

627 

628 

629 

630 

631 

632 GO 59 

633 

634 

635 

636 

637 

638 1 

639 1 

640 1 

641 GO 60 

642 

643 GO 61 

644 

645 

646 

647 

648 

649 
630 

651 1 

652 1 

653 1 

654 GO 62 

655 

656 

657 

658 

659 

660 
661 
662 

663 

664 

665 

666 
667 
66 8 1 

669 1 

670 1 


FMP E0037 A 

STU 9054 

RAU H0036 B 

NZU GO 59 

FMP F0038 A 

FMP R0038 A GO 59 

FSB 9054 

FDV M0032 C 

STU F0038 A 

NZA GO 60 

AXC 0001 

AXA 0001 GO 58 

FLUX CALC POINTS 37-32 

AXA 0005 

SXC 0006 GO 61 

FMP M0077 C 

STU 9000 

RAU F0032 A 

FSB 9000 

STU F0032 A 

NZA GO 62 

SXC 0001 

SX A 0001 GO 61 

FLUX CALC POINT 31 

AXC 0005 

RAU F0033 

FMP M0040 C 

STU 9000 

RAU F0031 

FSB 9000 

STU F 0031 

RAU F0032 

FMP M0071 C 

STU 9000 

RAU F0031 

FSB 9000 

STU F 0031 

AXA 0008 GO 63 

FLUX CALC POINTS 30-22 


671 

GO 63 FMP 

M0070 

c 



672 

STU 

9000 




673 

RAU 

F 0022 

A 



674 

FSB 

9 000 




675 

STU 

F0022 

A 



676 

NZA 



GO 

64 

677 

SXC 

0001 




678 

SX A 

0001 


GO 

63 


0650 39 2834 0234 
0234 21 9054 0441 
0441 60 4752 0207 
0207 44 0261 0162 
0261 39 3766 0266 
0266 39 3647 0162 
0162 33 9054 0491 
0491 34 6867 0167 
0167 21 3766 0269 
0269 40 #072 0323 
0072 58 0001 0128 
0128 50 0001 0650 


0323 50 0005 0229 
0229 59 0006 0235 
0235 39 6912 0212 
0212 21 9000 0319 
0319 60 3760 0365 
0365 33 9000 0545 
0545 21 3760 0313 
0313 40 0316 0217 
0316 59 0001 0122 
0122 51 0001 0235 


0217 58 0005 0523 
0523 60 1761 0615 
0615 39 6875 0275 
0275 21 9000 0583 
0583 60 1759 0663 
0663 33 9000 0443 
0443 21 1759 0312 
0312 60 1760 0665 
0665 39 6906 0506 
0506 21 9000 0314 
0314 60 1759 0364 
0364 33 9000 0493 
0493 21 1759 0362 
0362 50 0008 0550 


0550 39 6905 0205 
0205 21 9000 0263 
0263 60 3750 0255 

0255 33 9000 0185 
0185 21 3750 0353 
0353 40 0256 0157 

0256 59 0001 0062 
0062 51 0001 0550 


CTT 


CY-5 back E-113 


35 



0157 58 0008 0414 
0414 60 6874 0379 
0379 39 1751 0402 
0402 21 9000 0309 
0309 60 1749 0653 
0653 33 9000 0633 
0633 21 1749 0452 
0452 60 1750 0455 
0455 39 6896 0346 
0346 21 9000 0104 
0104 60 1749 0154 
0154 33 9000 0384 
0384 21 1749 0502 
0502 50 0008 0600 


0600 39 6895 0395 
0395 21 9000 0403 
0403 60 3740 0445 
0445 33 9000 0225 
0225 21 3740 0343 
0343 40 0146 0247 
0146 59 0001 0252 
0252 51 0001 0600 


0247 58 0008 0204 
0204 60 6917 0421 
0421 39 1741 0591 
0591 21 9000 0299 
0299 60 1739 0543 
0543 33 9000 0573 
0573 21 1739 0192 
0192 60 1740 0396 
0396 39 6886 0286 
0286 21 9000 0593 
0593 60 1739 0643 
0643 33 9000 0623 
0623 21 1739 0242 
0242 50 0009 0100 


0100 39 6885 0035 
0035 21 9000 0043 
0043 60 3729 0083 
0083 33 9000 0013 
0013 21 3729 0032 



36 


732 




NZ A 


GO 

68 



0032 

40 

0085 

0036 

733 




sxc 

0001 





0085 

59 

0001 

0291 

734 




SX A 

0001 

GO 

67 



0291 

51 

0001 

0100 

H5 


GO 

6 8 

AXC 

0009 





0036 

58 

0009 

0042 

736 




RAU 

<CHEK 





0042 

60 

9018 

0341 

737 




BM I 

GO 88 





0341 

46 

0244 

0245 

738 




RAL 

CODE 





0245 

65 

0198 

0203 

739 




LDD 

GO 88 

PUNCH 



0203 

69 

0244 

1930 

740 


GO 

88 

NZC 

GO 80 

GO 

79 



0244 

48 

0050 

0248 

741 

1 











742 

1 




STORE 

+ PUNCH NORMAL FLUX 






743 

1 













744 


GO 

79 

RAU 

KCHEK 





0248 

60 

9018 

0630 

745 




3MI 

GO 80 





0630 

46 

0050 

0580 

746 




RAU 

F 000 1 





0580 

60 

1729 

0440 

7 4? 




STU 

NORMF 





0440 

21 

0074 

0490 

748 




RAL 

PNORM 





0490 

65 

0680 

0540 

7 49 




LDD 

GO 80 

PUNCH 



0540 

69 

0050 

1930 

7 50 

1 











751 

1 




NEUTRON 

REGENERATION 






752 

1 













753 


GO 

80 

SX A 

0010 

GO 

69 

REGTON 

I 

0050 

51 

0010 

0006 

7 54 


GO 

69 

RAU 

V 000 9 

B 




0006 

60 

5692 

0047 

7 5 5 




NZU 


GO 

71 



0047 

44 

0001 

0002 

7 3 6 




FMP 

F0011 

A 




0001 

39 

3739 

0039 

73? 




FAD 

booh 

A 




0039 

32 

3777 

0003 

7 5 8 




STU 

booh 

A 




0003 

21 

3777 

0030 

7 5 




NZ A 


GO 

81 



0030 

40 

0033 

0034 

7 6 C 


GO 

81 

AX A 

0001 

GO 

69 



0033 

50 

0001 

0006 

7 6 1 


RAU 

KCHEK 





0034 

60 

9018 

0041 

7 6 2 




BM I 

GO 82 





0041 

46 

0044 

0045 

7 6 3 




RAL 

PCB1 





0045 

65 

0048 

0053 

764 




LDD 

GO 82 

PUNCH 



0053 

69 

0044 

1930 

765 

1 












766 


GO 

82 

SX A 

0010 

GO 

71 

REGION 

II 

0044 

51 

0010 

0002 

76 7 


GO 

71 

RAU 

V0018 

B 




0002 

60 

5701 

0005 

76 8 




NZU 


GO 

73 



0005 

44 

0009 

0010 

769 




FMP 

F0021 

A 




0009 

39 

3749 

0049 

7 7 r 




FAD 

B 002 2 

A 




0049 

32 

3788 

0015 

771 




STU 

BOO 22 

A 




0015 

21 

3788 

0091 

7 72 




NZ A 


GO 

83 



0091 

40 

0094 

0095 

773 




AX A 

0001 

GO 

71 



0094 

50 

0001 

0002 

"*74 


GO 

83 

RAU 

KCHEK 





0095 

60 

9018 

0141 

7 7 5 




BM I 

GO 84 





0141 

46 

0144 

0145 

776 




RAL 

PCB2 





0145 

65 

0098 

0103 

777 




LDD 

GO 84 

PUNCH 



0103 

69 

0144 

1930 

778 

1 












77c 


GO 

84 

SX A 

0010 

GO 

73 

REGION 

III 

0144 

51 

0010 

0010 

7 8 C 


GO 

73 

RAU 

V002 7 

B 




0010 

60 

5710 

0065 

78 1 




NZU 


GO 

93 



0065 

44 

0019 

0020 

782 




FMP 

F 00 3 1 

A 




0019 

39 

3759 

0059 

7 83 




FAD 

B 0033 

A 




0059 

32 

3799 

0025 

7 86 




STU 

B 003 3 

A 




0025 

21 

3799 

0052 


E-113 


E-113 


37 


785 



NZ A 


GO 85 


0052 

40 

0055 

0056 

786 



AXA 

0001 

GO 73 


0055 

50 

0001 

0010 

787 


GO 85 

RAU 

KCHEK 



0056 

60 

9018 

0191 

788 



BMI 

GO 86 



0191 

46 

0194 

0195 

789 



RAL 

PCB3 



0195 

65 

0148 

0153 

790 



LDD 

GO 75 

PUNCH 


0153 

69 

0106 

1930 

791 


GO 93 

RAU 

KCHEK 



0020 

60 

*018 

0241 

792 



BMI 

GO 86 

GO 75 


0241 

46 

0194 

0106 

793 

1 










794 

1 


NORMALIZE FLUX 






795 

1 










796 


GO 75 

ral 

CODE 



0106 

65 

0198 

0453 

797 



LDD 


PUNCH 


0453 

69 

0306 

1930 

798 



RSA 

0037 

GO 76 


0306 

81 

0037 

0112 

7 99 


GO 76 

RAU 

F0038 

A 


0112 

60 

3766 

0171 

POO 



FDV 

NORMF 



0171 

34 

0074 

0124 

801 



STU 

N 0038 

A 


0124 

21 

3839 

0092 

802 



NZ A 


GO 77 


0092 

40 

0495 

0196 

803 



AXA 

0001 

GO 76 


0495 

50 

0001 

0112 

804 


GO 77 

RAL 

PCNF 



0196 

65 

0249 

0503 

805 



LDD 

GO 86 

PUNCH 


0503 

69 

0194 

1930 

806 

1 










807 

1 


GO TO 1 

NEXT GROUP OR 

CALC POWER 





808 

1 










809 


GO 86 

NZB 


MOVEB 


0194 

42 

0397 

0448 

810 



AXB 

0001 



0397 

52 

0001 

0254 

811 



AXC 

0082 

GO 78 


0254 

58 

0082 

0463 

812 

1 










813 

1 




CONSTANTS 






814 

1 










815 


9001 

10 

0000 

0051 

ENERGY 

9001 

10 

0000 

0051 

816 


9002 

20 

OOOO 

0051 

GROUP CODE 

9002 

20 

0000 

0051 

817 


9003 

30 

0000 

0051 


9003 

30 

0000 

0051 

818 


9004 

40 

0000 

0051 


9004 

40 

0000 

0051 

819 


9005 

50 

0000 

0051 


9005 

50 

0000 

0051 

820 


9006 

60 

0000 

0051 


9006 

60 

0000 

0051 

821 


9007 

70 

0000 

0051 


9007 

70 

0000 

0051 

822 


9008 

80 

0000 

0051 


9008 

80 

0000 

0051 

823 


9009 

90 

0000 

0051 


9009 

90 

0000 

0051 

824 

1 










825 


MIN 1 

- 10 

0000 

0051 

OTHER 

0466 

-10 

0000 

0051 

826 


CT 1 2 

12 

0000 

0052 

CONSTANTS 

0579 

12 

0000 

0052 

827 


CT 38 

38 

0000 

0052 


0280 

38 

0000 

0052 

828 


CT 39 

39 

0000 

0052 


0388 

39 

0000 

0052 

829 


CT40 

40 

0000 

0052 


0679 

40 

0000 

0052 

830 


CT 96 

96 

0000 

0052 


0542 

96 

6000 

0052 

831 

1 










832 


Q0001 

00 

000 

000 

SIGMA 0 

1574 

00 

0000 

0000 

833 


Q0010 

00 

000 

000 

FOR FIRST 

1583 

00 

0000 

0000 

8 34 


Q0019 

00 

000 

000 

GROUP 

1592 

00 

0000 

0000 

835 


00028 

00 

000 

000 


1601 

00 

0000 

0000 

836 

1 










837 


CODE 

06 

9409 

0001 

PUNCH 

0198 

06 

9409 

0001 



38 


838 


PCB1 

06 

B0001 

0011 

CONSTANTS 

0048 

06 

1767 

0011 

839 


PCB2 

06 

B0012 

0011 



0098 

06 

1778 

0011 

840 


PCB3 

06 

B0023 

0011 



0148 

06 

1789 

0011 

841 


PCNF 

06 

N 000 1 

0038 



0249 

06 

1802 

0038 

842 


PCR 

06 

R 0001 

0038 



0581 

06 

1610 

0038 

843 


PCKEF 

06 

9012 

0002 



0662 

06 

9012 

0002 

844 


PN0RM 

06 

NORMF 

0001 



0680 

06 

0074 

0001 

845 

1 











846 


C0082 - 

40 

0000 

0051 

01 05 

09 

1801 

-40 

0000 

0051 

847 


E0038 

00 

000 

000 

E38 


0835 

00 

0000 

0000 

848 

1 











849 


DEL K 

10 

0000 

0047 

DEL K 

TEST 

0625 

10 

0000 

0047 
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TABLE I. - STORAGE LOCATIONS OF INPUT 


Storage 

region 

Storage 

location 

Input value 

A 

0753 to 0761 

a a> a 'b> •** a l 

D 

0762 to 0770 

• D I,a> D I,b> **• D I,i 


0771 to 0779 

D II,a' D II/b> **• D II,i 


0780 to 0788 

D III,a> D III,b> D III,i 


0789 to 0797 

D IV,a-> VrVj'b-' •** D IV,i 

Q 

1575 to 1582 

Z q,I,a> Z q,I,b> ••• Z q,I,h 


1584 to 1591 

Z q,II,a' Z q,II,b> ••* Z q,II,h 


1593 to 1600 

2 q,III,a' Z q,III,b' •** Z q,III,h 


1602 to 1609 

Z q,IV,a' Z q,IV,b> ••‘ Z q,IV,h 

s 

1648 to 1656 

•** ^i 


1657 to 1665 

Z II,a> Z II,b> •** Z II,i 


1666 to 1674 

^III,a' Z III,V *** Z IH,i 


1675 to 1683 

Z IV,a-> Z IV,b-> •** Z IV,i 

V 

! 1684 to 1692 

VZ f,I,a' V f,I,b' V2 f,l,i 


1693 to 1701 

v2 f,II,a' vZ f,II,b' vZ f,II,i 


1702 to 1710 

vZ f,III,a> v2 f,III,b> ••• v2 f,III,i 

None 

I 

9014 to 9017 

h I ,h II ,h III ,h IV 
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TABLE II. - ORDER OF PUNCHED OUTPUT 


Program 

Output value 

Location 

Remarks 

instruc- 


punched 


tion 


out 


489,490 

k e f f ,Ak e ff 

9012,9013 

(Punched for every iteration) 

739,740 

Group code number 

(9009+IAB)* 


748,749 

Normal flux. cp_, 

l,a 

i 

0680 


763,764 

y ' ^I,n,k 

1767 to 

For points 1 to 11, when 


n=a 

1777 

vZ f,I,n f 0 

776,777 

^ ^II,n,k 

1778 to 

For points 11 to 21, when 


n=a 

"J 

1788 

vZ f,II,n £ 0 

789,790 

/ . P III,n,k 

1789 to 

For points 21 to 31, when 


n==a 

1799 

v2 f,III,n i 0 

796,797 

Group code number 

(9009+IAB)* 


804,805 

^l,n ^38, n 

1802 to 


*L,a 

1839 


494,495 

r l> • • * r 38 

1610 to 



164 7 


496,497 

k ef f,Ak e ff 

9012,9013 

Last values punched out 


" When the first neutron energy group is being considered in the 
program, the value of IAB equals -0008; this results in the 
punching of the contents of location 9001, which equals 
10 0000 0051. Here, the last two digits indicate the position 
of the decimal to be 1.0 (ref. 4). 
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